Hydrogen site occupancy and strength of forces in nano-sized metal hydrides 
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The dipole force components in nano-sized metal hydrides are quantitatively determined with 
curvature and x-ray diffraction measurements. Ab-initio density functional theory is used to cal- 
culate the dipole components and the symmetry of the strain field. The hydrogen occupancy in a 
100 nm thick V film is shown to be tetrahedral with a slight asymmetry at low concentration and a 
transition to octahedral occupancy is shown to take place at around 0.07 [H/V] at 360 K. When the 
thickness of the V layer is reduced to 3 nm and biaxially strained, in a Feo.sVo.s/V superlattice, the 
hydrogen unequivocally occupies octahedral 2-like sites, even at and below concentrations of 0.02 
[H/V]. 



PACS numbers: 68.55. Ln, 61.05.cm, 61.05.cp, 68.65. Cd 
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I. INTRODUCTION 

Hydrogen and its isotopes can be absorbed in large 
quantities into the interstitial sites of many transition 
metals (iii Hydrogen constitutes the simplest impurity 
problem and exhibits many properties that can be ex- 
plained by a lattice gas modelj^i^ The interaction be- 
tween the absorbed hydrogen atoms is of special inter- 
est, since it is mediated by the host metal through elas- 
tic deformations <^i^ For example, the long ranged attrac- 
tive interaction is responsible for the condensation into 
liquid- and solid like phases. Understanding and quan- 
tifying the forces exerted by hydrogen is thus crucial for 
the understanding of phase transitions in metal hydrides. 
Due to the long range of the inter-hydrogen attraction in 
e.g. the transition metals, the phase formations addi- 
tionally depend crucially on the elastic boundary condi- 
tions. Thus when studying hydrogen in nano-sized sys- 
tems, large changes are expected when the aspect ratio 
and shape are altered. 

Measurements of the hydrogen induced local elastic 
stress have been carried out in bulk metal systemsi^i 
by measuring the components of the dipole force tensor, 
which is defined by 
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where Ri is the i-th component of the position of the 
77i-th metal atom and /'" is the j-th component of the 
force exerted by hydrogen on metal atom m. When hy- 
drogen atoms (concentration c in atomic ratio [H/M]) are 
distributed randomly in the host structure, the average 
stress atj can be related to the components of the dipole 
tensor and the atomic volume of the metal (fi). The 
dipole force tensor is the connecting quantity between 



the atomistic and the elastic-continuum description of 
lattice distortions originating froin point defects in crys- 
tals.— The dipole force tensor is analogous to an electric 
or magnetic dipole and has been used for instance in the 
description of ferroelasticitj«S. 

In bee metals two distinct high-symmetry sites can be 
occupied by hydrogen, tetrahedral (T) and octahedral 
(O). These sites exhibit tetragonal symmetry and the 
dipole tensor can be written for the case of Oz or T^ sites 
as 
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The dipole force tensor for other hydrogen absorption 
sites are straightforward permutations of Eq. [2] In gen- 
eral, the components of the dipole tensor are different 
for hydrogen occupying the two sites and thus serve as a 
signature for the type of site occupied. 

The hydrogen-induced stress field in nano systems has 
previously been investigated by KirchheintiSrJ^ and co- 
workers, where they showed, for instance, that the clastic 
constants of a Nb film could be determined from mea- 
surements of stress and hydrogen concentration. Stress 
measurements have also been used to follow the forma- 
tion of different hydride phases in polycrystalline films of 
YHjjiii Hydrogen absorbed into low-dimensional struc- 
tures such as epitaxial thin films and superlattices can ex- 
perience different values of the dipole components due to 
strain effects. For example Fe/V superlattices are tetrag- 
onally strained with a ratio that depends on the ratio of 
the constituent materials.—""— The hydrogen absorption 
properties of thin films and superlattices, investigated in 
the past, have shown distinctly different properties than 
that of the bulkii&""— However, it has never been convinc- 
ingly demonstrated into what site hydrogen goes in these 



types of structures nor to what extent the finitencss of 
the extension dictates the evolution of the stress field. 

We measure directly the dipole components in a vana- 
dium film and a Feo.s/Vo.s (001) superlattice, as a func- 
tion of the hydrogen concentration and determine the 
site occupancy in these types of structures. We will 
also determine the components of the dipole force ten- 
sor from first principles by calculating the Hellmann- 
Feynnian forces using density functional theory. This 
allows us to verify the occupancy of the two distinct sites 
as a function of H concentration, which were experimen- 
tally determined from measurements of the dipole com- 
ponents. 
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II. EXPERIMENTAL DETAILS 
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A. Sample preparation and characterization 



The Feo.sVo.s/V (6/21) ML superlattice with 100 rep- 
etitions was grown on a 20x20x1 mm^ polished single- 
crystal MgO (001) substrate using UHV-based mag- 
netron co-sputtering. The base pressure in the system 
was 2xl0~^ Pa prior to deposition with the dominating 
residual gas being hydrogen. The sample was deposited 
from targets of iron and vanadium onto the substrate 
held at 573 K. Argon, with purity better than 99.9999 
% at a pressure of 0.27 Pa was used as the sputtering 
gas. The deposition rates of Fe and V were 0.0134 nm 
s~^ and 0.0143 nm s"""^, respectively. The first vanadium 
layer was directly deposited onto the MgO substrate and 
the superlattice was capped with a vanadium layer fol- 
lowed by a 10 nm palladium layer. The palladium was 
sputtered from a target with 99.99 % purity after cooling 
the sample to room temperature. The palladium layer 
acts as a catalyst for the hydrogen dissociation and pro- 
tects the underlying structure from oxidation. The 100 
nm vanadium film was grown epitaxially on MgO (001) 
10x10x0.5 mm substrate using the same recipe as above 
except a 5 nm palladium cap was grown on top of the 
film. 

The 100 nm V film was measured with x-ray diffraction 
and has a coherence length of 50 nm in the z direction and 
a mosaicity of 0.2 degrees. Thus the film can be viewed 
as a single-crystal since the coherence length is of the 
order of the film thickness. The Fco.sVq.s/V superlattice 
structure has been characterized previously using x-ray 
and neutron refiectivity as well as x-ray diffraction and 
transmission electron microscopy i^ It was determined to 
be a single-crystal that showed no signs of changes in 
the unloaded XRD or reflectivity pattern during repeated 
cycles of hydrogenation. This is interpreted as an absence 
of subsequent structural degradation. 



FIG. 1: A schematic illustration of the curvature measure- 
ment. A laser beam is split into a rectangular array of beams 
that reflect from a curved sample and is detected by a CCD 
camera. The change in spot spacing can be related to a change 
in curvature. 
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FIG. 2: A typical array of spots reflected from the sample 
surface and detected by the camera. The spot spacing yields 
the curvature according to Eq. (3] 



B. Curvature measurements 

A multi-beam optical stress sensor (MOSS) is used to 
measure curvature in real-time, which is schematically 
depicted in Fig. [T] A laser beam is split into multiple 
beams using two etalons, which reflect from the sample 
surface into a CCD camera. The whole laser system is 
enclosed and kept at a constant temperature to mini- 
mize the influence of thermal drift on the measurements. 
Commercial software (kSpace associates) is used to an- 
alyze the reflected pattern and calculate the curvature. 



The curvature instrument includes an ultra-high vacuum 
compatible chamber, which can be heated to 400 K. Ca- 
pacitive membrane gauges are used to measure the pres- 
sure in the range 0.01-10 Pa. Figure [5] shows a typical 
image recorded by the CCD camera from the sample. 
The change in spot spacing d is directly related to the 
curvature k by 
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where L is the distance from the laser to the sample and 
a is the angle of incidence (see Fig. [1]). The curvature of 
the film is measured directly as a function of temperature 
and external H2 pressure. The curvature of the substrate 
is related to the stress in the film by Stoney's formula 
adopted for substrates with a (001) orientation: 



where ei is shorthand notation for en (Voigt notation), 
etc., Sij are the compliance constants of the crystal and 
have units of inverse of pressure. The compliance con- 
stants for vanadium were converted from [Ref. |24| and 
are sn = 6.821 x IQ-^^ p^-i and S12 = -2.336 x IQ-^^ 
Pa^^. Any change in the compliance constants due to 
hydrogen absorption and temperature can be completely 
neglected in the concentration and temperature range in- 
vestigatedi^ However, since the superlattice is tetrago- 
nally strained and contains V, which is affected by the 
proximity to Fe, the elastic constants will be slightly dif- 
ferent compared to the pure V constants. This difference 
is neglected in the present study. 

The strain in the different directions for a free sample 
can now be evaluated by inserting Eq. ^ and Eq. ([2]) 
into Eq. (O, which gives: 
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ei = u[si2A+{sii + si2)B] 

£2 == 7i[si2A+{sii + Si2)B] 

£3= f^[snA + 2si2B]. 



(6) 



where s,, are the compliance constants of the MgO sub- rj^j^^ ^^j^^^ ^^^^^^ -^ ^-^^^ evaluated as 



strate and <sub and ifiim are the thickness of the substrate 
and the film, respectively. A didactic derivation of this 
version of Stoney's formula can be found in Ref. [23| . 
The stress calculated from Stoney's formula is positive 
when the substrate curves as shown in the right panel of 
Fig. \E\ (concave).— The resolution of the curvature mea- 
surement is approximately 2 x 10~^ m~^ and a typical 
uncertainty in stress is of the order of 0.04 GPa. 



C. In situ X-ray diffraction 

A specially designed UHV scattering chamber was used 
for determining both the hydrogen induced expansion 
and structural changes in the samples. This scattering 
chamber allows in-situ exposure to hydrogen in a wide 
temperature and pressure range, and was mounted on a 
Bruker Discover D8 X-ray diffractometer equipped with a 
parallel X-ray beam (CuKai A=0. 15406 nm). The sam- 
ple was exposed to highly purified H2-gas at different 
pressures at constant temperature (0.01-10 Pa). Thermo- 
dynamic equilibrium was ensured by continuously mon- 
itoring the (002) Bragg peak for changes after a sudden 
pressure increase. 



D. Dipole components of a clamped film 

The strain propagating through the metal due to the 
dipole force tensor is given by elasticity theory as 
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To understand what happens to the volume change when 
a film is clamped on a substrate, the following thought 
experiment can be invoked. Assume that the film fits 
perfectly on the substrate and is free to expand according 
to Eq. ([6|). Next, we impose an imaginary external stress 
in the plane of the film, which completely negates the in- 
plane expansion. We do not know what stress is required 
but the strain has to be 
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By using Eq. ([5]), we can immediately write down some 
general relations for this kind of application of stress: 
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Combining the relations in Eq. (jlOp yields 
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The total volume change due to hydrogen in a clamped 
film is thus 
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E. Concentration conversion 

Sicvert's law was used to convert the H2 pressure p and 
temperature to concentration 



c = K'\T). 



(14) 



where K is the Sievert constant. For the 100 nm V 
fihn the constant was obtained from the hterature on 
50 and 10 nm fihns, which are consistent with bulk solu- 
bility^!^ 

The Sievert constant for the Feo.sVo.s/V superlattice 
was obtained by measuring the relative change in the re- 
sistivity as a function of the applied hydrogen pressure 
at different temperatures. The relation between resistiv- 
ity changes and concentration was obtained from previ- 
ous neutron reflectivity measurements j^ The value of the 



Sievert constant was determined to be K^ 
at 360 K. 
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FIG. 3: Illustration of the lattice distortion induced by 
a hydrogen atom (small dark ball) at the interstitial tetra- 
hedral site in bcc bulk vanadium (larger grey balls). The 
overlayer shows the magnitude of the displacement field as 
changes in the inter-vanadium bond length (white: contrac- 
tion; black: expansion). The contour lines of constant dis- 
tortion are shown as white lines. The far field expansion is 
isotropic with H at the tetrahedral site. 



Theoretical methods 



The scalar-relativistic ab-initio density functional the- 
ory (DFT) calculations were performed using the projec- 
tor augmented wave (PAW)21i2^ method as implemented 
The exchange-correlation interaction was 
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treated in the generalized gradient approximation (GGA) 



in the parameterization of Perdcw, Burke, and Ernzer- 
hof (PBE)i^ A cut-off energy of 500 cV was used and a 
Gaussian smearing with a width of cr = 0.05 eV for the 
occupation of the electronic levels. For the bulk calcu- 
lations the 3p^, 3d^, and 4^^ states of vanadium were 
treated as valence states. In the superlattice calculation 
the Sd'^ and 4s^ states of vanadium as well as the Sd"^ 
and 4s^ states of iron were treated as valence states. 

The bulk bcc vanadium was modeled by a 3 x 3 x 3 (54 
atoms) or a 4 X 4 X 4 (128 atoms) supercell, respectively, 
with a calculated lattice constant of 3.007 A. The iron - 
vanadium superlattice was modeled by repeating a 2 x 2 
supercell, consisting of 7 monolayers (ML) of iron and 7 
ML of vanadium, in the z-direction. In the experiments 
the superlattice was grown epitaxially on a magnesium 
oxide (001) surface. This compresses the in-plane lattice 
constant of vanadium by 1.61%. In the calculations, we 
simulated this relative change by starting from the cal- 
culated lattice constant of MgO (4.21 A), which leads to 
an in-plane lattice constant of V of 2.93 A. Note: the cal- 
culated bulk lattice constant of V with the smaller set of 
valence states is 2.98 A. Since the alloy does not take up 
appreciable hydrogen at the temperatures and pressure 
used in this study, it was deemed sufficient to use pure 
Fe instead. The full size of the superlattice was too large 
to be modeled within a reasonable time frame so Fe/V 
7/7 ML was chosen as a reasonable compromise between 
ease of conrputation and still being physically relevant. 

A Monkhorst-Pack F-centered 4x4x4 k-point mesh 
(36 k-points in the irreducible wedge of the Brillouin- 
Zone) was used for the structural relaxations of the 3 x 
3x3 bcc V supercell, a 3 x 3 x 3 k-point mesh for the 
4x4x4 bcc V supercell, and a 7 x 7 x 3 k-point mesh 
for the Fe/V superlattice. Spin-polarization was taken 
into account in the bulk calculations. The relaxation 
cycle was stopped when the Hellmann-Feynman forces 
had become smaller than 5 ■ 10~^ eV/A. 



III. RESULTS AND DISCUSSION 

A. Stress and strain 

Figure [3] shows calculated local strain fields in bulk bcc 
vanadium with hydrogen occupying a tetrahedral z site 
(T2). The overlayer shows the magnitude of the displace- 
ment field as changes in the inter- vanadium bond length 
(white: contraction; black: expansion). Within the nu- 
merical accuracy, the displacement field is isotropic. Fig- 
ure [3] shows the strain field with hydrogen occupying a 
octahedral z site (O^). Notice the strong asymmetry 
compared to the tetrahedral case. The forces exerted on 
the vanadium lattice before the lattice is allowed to relax 
can be used to approximate the dipole components ac- 
cording to Eq. [TJ The results are shown in Tab. U along 
with experimental values from the literature and experi- 
mental values from curvature and x-ray diffraction mea- 
surements. It is important to note that when extract- 
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FIG. 4: Illustration of the lattice distortion induced by a hy- 
drogen atom (small dark ball) at the interstitial octahedral 
site in bcc bulk vanadium (larger grey balls). The overlayer 
shows the magnitude of the displacement field as changes in 
the inter-vanadium bond length (white: contraction; black: 
expansion). Additionally, in white the contour lines of con- 
stant distortion are shown. Note that the expansion has a 
rotation symmetry with respect to an axis parallel to the z- 
axis that passes through H atom. 
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FIG. 6: In-plane stress as a function of concentration. The 
solid points correspond to the Feo.sVo.s/V superlattice and 
the open points represent the 100 nm V film. The data is 
taken at 360 K. 




FIG. 5: Curvature changes due to film stress. The red arrows 
signify the stress experienced by the substrate due to the film. 
Panel (a) illustrates how compressive stress leads to convex 
curvature. Panel (b) shows that tensile in-plane stress results 
in concave curvature. 



ing dipole components from experimental data, the elas- 
tic constants are assumed not to be affected. In other 
words, changes in the elastic compliance of the mate- 
rial due to hydrogen is incorporated into the values of 
the dipole components. Hovifcver, the calculations yield, 
within the approximations of DFT, estimates of the phys- 
ical, Hellmann-Feynman, forces. 

A large discrepancy between calculated and measured 
dipole components from literature is found. To verify the 
computational approach, the calculations were repeated 
for one hydrogen atom in niobium and the results were in 
good agreement with previous calculations as well as ex- 



perimental data. Furthermore, in the present calculation 
of H absorbed in bulk vanadium, a larger supercell, i.e. 
containing 128 atoms, and a larger set of valence states 
have been employed than in previous studies. Hence, our 
results are well-converged with respect to these param- 
eters and the problem might arise from the vanadium 
pseudo-potentials as such. This conclusion is supported 
by a recent theoretical investigation^, which also failed 
to correctly predict the elastic constants of vanadium at 
the GGA level. Even though the calculations fail to pre- 
dict quantitatively the values of the dipole components, 
they qualitatively predict the correct symmetry of the 
strain field for the two different sites. The calculations 
show a decrease in the lattice spacing in the x-y directions 
when hydrogen occupies the Oz site as shown in Fig. UJ in 
agreement with experimental diffraction experiments)^ 
We are now in a position to describe how hydrogen 
absorption into a film grown on a substrate dictates the 
curvature of the substrate. For a quantitative deriva- 
tion see section IIIDI If hydrogen occupies the tetrahe- 
dral sites, the strain field propagates to the boundaries 
and contributes to an isotropic expansion. The substrate 
acts as a strong restoring force on the in-plane expan- 
sion, which leads to a compressive stress at the surface 
of the substrate. This in turn produces a bending mo- 
ment and the substrate bends convex as shown in panel 
(a) of Fig. [5] Therefore, it is expected a priori that a 
fiat sample, where hydrogen occupies tetrahedral sites, 
exhibits a convex curvature. When hydrogen occupies 
octahedral z sites in the bulk, as found in the /3-phase of 



VHa;, the material in the bulk contracts in the plane and 
expands out of plane. To fit the sample on a substrate 
the sample has to be stretched in the plane with a tensile 
stress. The restoring force of the substrate is thus tensile 
and the sample bends concave as shown in panel (b) of 
Fig. \5\ Therefore it is expected that a flat sample with 
hydrogen occupying octahedral z sites exhibits a concave 
curvature. 
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FIG. 7: Changes in the total energy of a Fe/V superlattice as 
a function of hydrogen concentration. The hydrogen atoms 
are absorbed at distorted tetrahedral (open squares) and oc- 
taliedral (full circles) sites. 



Figure El shows the in-plane stress of the 100 nm vana- 
dium film (open circles) as a function of concentration. 
The stress is obtained by converting the measured cur- 
vature using Eq. 21 A sharp change in the stress is ob- 
served at 0.065 in [H/V]. Below this transition the stress 
increases linearly with concentration as expected from 
the definition Eq. ([1]). From the discussion above, we 
therefore assign the region below the transition to hy- 
drogen occupying sites that have isotropic dipole compo- 
nents (tetrahedral sites) . The region above the transition 
is where the dipole components are sufficiently different 
from each other to cause a contraction in the plane of 
the sample. Thus from a qualitative point of view, these 
results are consistent with tetrahedral occupancy at low 
concentrations and octahedral at concentrations above 
0.065 [H/V], which resembles the bulk occupancies. 

Fig.[n]also shows the stress for a Fco.sVo.s/V superlat- 
tice (filled circles). In stark contrast to the film results, 
we see an almost immediate increase (positive) in the 
stress with increasing concentration. There are therefore 
significant differences in the response from a film and 
a superlattice, which will be discussed more thoroughly 
below. 

Using the dipole components from literature, Eq. [T3] 
predicts the stress on a MgO substrate to be ai/c = 
4.9904 GPa for octahedral z sites and cti/c = -14.0913 
GPa for tetrahedral sites. We see that the two stresses 



have different signs, which is expected from the net con- 
traction in the plane for hydrogen occupying octahedral 
sites and a net expansion for hydrogen occupying tetra- 
hedral sites (see also Fig. [5|). 

Table H] shows literature values for the sum and the dif- 
ference of the dipole components for tetrahedral sites T 
and Oz sites. The volume change of a free sample can be 
calculated using Eq. [T] The volume change expected from 
a clamped film is also calculated from Eq. [121 Notice that 
the volume change is larger for the clamped film when 
hydrogen occupies Oz sites than its free counterpart. In 
a previous publication^^, we used a simplified version of 
the above equation, which yielded a lower than expected 
value based on octahedral occupancy. Using the full elas- 
tic theory the results previously published can now be 
understood and are consistent with octahedral z site oc- 
cupancy. Note that the out of plane expansion does not 
depend on the in-plane components B. However the in- 
planc stress measured by the curvature change depends 
on both A and B. By measuring the out-of-plane expan- 
sion the A components can be uniquely determined, as 
described in the following chapter. 



B. Lattice expansion 

Figure \E\ shows the relative change in the doo2 lattice 
distance as a function of pressure (upper panel) and con- 
centration (lower panel) at 400 K for the 100 nm vana- 
dium film. The lower panel shows the data converted 
from pressure to concentration in the region below the 
transition. The slope is 0.1189(7), which gives, using 
Eq. [T^ a dipole component A = 2.35(1) eV. Using the 
value of A obtained from the expansion and Eq. [T3] we 
can determine B by comparing to Fig. Ej This is shown 
by the solid line in the concentration range 0.00-0.068 
[H/V]. The component is slightly smaller than A with 
a value B = 2.0(1) eV. The horizontal line in the up- 
per panel of Fig. [8] gives using Eq. [12] with c = 0.44 a 
value A — 3.8(1) cV, which is consistent with octahe- 
dral occupancy. The concentration corresponds to the 
end of the isotherm plateau at 400 K. It is worth noting 
that large changes in pressure are necessary to change 
the concentration beyond the plateau, which is why a 
single concentration was used. For the region after the 
transition Eq. [T2] was used assuming an asymmetry of 
the components with A = 3.8(1) eV as determined from 
diffraction and B = 1.7(2) eV. 

It can thus be unambiguously concluded that the hy- 
drogen occupies tetrahedral sites in 100 nm thick V films 
below the phase boundary and octahedral sites beyond 
the co-existence region. This film can be viewed as a con- 
sistency check for the methodology presented. The dipole 
components are in quantitative agreement with literature 
except for a slight tetragonality in the B component for 
tetrahedral occupancy. This implies that Snoek anelastic 
relaxatioii^ could in principle be found for hydrogen in 
vanadium films, which has previously been excluded for 
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0.21(2) 


Palsson (2010) 


Fe/V 7/7 (calc) 


T 


12.97 


0.29 


5.15 


3.91 


- 


0.105 


Present work 


Fe/V 7/7 (calc) 


0. 


13.31 


1.95 


10.21 


1.55 


- 


0.140 


Present work 


Nb (calc) 


T 


10.45 


0.07 


3.31 


3.57 






Present work (Nb54Hl) 


Nb (calc) 


T 


11.0 


0.08 


3.47 


3.77 






Sundell (2004) 


Nb (calc) 


o. 


13.2 


2.11 


10.6 


1.3 






Sundell (2004) 



TABLE I: Dipole force tensor components for vanadium hydride. The components are calculated from the trace and the 
anisotropy. The expansion is calculated using A and B and Eq. ((T)). The clamped volume change is calculated from A using 
Eq. (|12|) . Items labeled as (calc) are results from ab-initio DFT calculations. 
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IV. CONCLUSIONS 



Since an unusually large in-plane stress was observed 
for the Feo.sVo.s/V superlattice, a fit to the data using 
Eq. [13] was not attempted. It is clear that although the 
sign from the stress measurements and the magnitude 
of the A component determined from previous neutron 
studies^ point towards an octahedral z site occupancy, 
we cannot explain the large stress observed. A calcu- 
lation of the change in total energy of a 7/7 Fe/V su- 
perlattice with hydrogen (Fig. [T]) shows that a distorted 
tetrahedral site is favored over the octahedral z site at 
low concentrations. It is possible that the occupancy of 
this distorted site is related to the unusual stress. Higher 
concentration measurements of the stress are required 
to elucidate this possibility which however awaits fur- 
ther study. The volume expansion in clamped vanadium 
films is larger than the free counter part, when hydrogen 
occupies the octahedral sites. This is due to the large 
asymmetry of the dipole components, as described by 
continuum elasticity theory. The present methodology 
is not restricted to hydrogen in metals but works for all 
types of interstitial defects that can be quantified by a 
dipole tensor. 



It is concluded that the force dipole components can 
be quantitatively measured by combining curvature and 
x-ray diffraction measurements in nano-systems. A 100 
nm film of vanadium loaded with hydrogen exhibits a 
stress field that is similar to tetrahedral occupancy be- 
low 0.065 [H/V] at 360 K. Above 0.068 [H/V] a sharp 
transition to octahedral sites occurs, which is consistent 
with the known phase boundaries of VHx- Hydrogen in 
superlattices composed of Fco.sVq.s/V are shown for the 
first time to occupy octahedral z-like sites at all measured 
concentrations. The symmetry of the strain fields due to 
hydrogen in different sites are qualitatively reproduced 
using ab-initio DFT calculations. The calculations also 
predict a change of site from tetrahedral to Oz at low 
concentrations. However, the quantitative discrepancies 
between our experiment and calculations show the need 
to treat the inter-atomic interactions in metal-hydrides 
with methods beyond the GGA level. The methodol- 
ogy presented here opens up new routes for investigating 
quantitatively, the evolution of local stress fields due to 
point defects in nano-sized materials. 
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